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Abstract 

I present improved predictions for the total hadronic cross section of stop-antistop production at 
hadron colliders including next-to-next-to-leading-order threshold corrections and approximated 
Coulomb corrections. The results are based on soft corrections, which are logarithmically en- 
hanced near threshold. I present analytic formulas for the NNLO scaling functions at threshold 
and explicit numbers for the total hadronic cross sections for the Tevatron and the LHC. Finally I 
discuss the systematic error, the scale uncertainty and the PDF error of the hadronic cross section. 
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Figure 1: LO production of a f, t* pair via gg annihilation (diagrams a-d) and qq annihilation (diagram e). 

1 Introduction 

The Minimal Supersymmetric Standard model (MSSM) is an attractive extension [TJ[2) of the 
very successful Standard Model. One property of this theory is its rich spectrum of new heavy 
particles which might be discovered at the LHC if they are lighter than sa 2TeV. Searches for 
superymmetric particles are performed at the Tevatron and the LHC. No superpartners have 
been discovered so far. In the mSugra framework of the Minimal Supersymmetric Standard 
Model (MSSM), the lightest stop squark, one of the scalar supersymmetric partner particle of the 
top quark, is assumed to be the lightest supersymmetric coloured particle, lighter than the other 
scalar quarks. This is due to large nondiagonal elements in the stop mixing matrix, see f. e. the 
review (3) . The lightest stop squark might be the first coloured SUSY particle to be discovered. 
The cross section delivers information about the stop mass or, if the mass of the stop squark 
is roughly known from elsewhere, information about its spin HJ. If these particles cannot be 
discovered at the Tevatron or LHC, precise cross sections help to improve mass exclusion limits. 
In this paper, I study the hadroproduction of stop-antistop-pairs 

VV l V p -> ht*X, i = l,2, (1) 

with its partonic subprocesses 

gg —> tjt* and qq —> iff, q = u,d,c,s,b, (2) 

including NNLO threshold contributions. The relevant leading order (LO) Feynman diagrams 
are shown in Fig. [TJ The production of mixed stop pairs ?if| or fat* starts at next to leading 
order (NLO) |5j and is therefore suppressed. This case will not be considered in this paper. 
The top parton density distribution in a proton is assumed to be zero, in contrast to the other 
quark parton density distributions. As a consequence, there is no gluino exchange diagram as 
for squark antisquark hadroproduction. For that reason, the qq channel is suppressed by a larger 

power of f> = yl — 4m? /s due to P-wave annihilation. The final state must be in a state with 
angular momentum I = 1 (denoted as P) to balance the spin of the gluon. Therefore, the case 
of stop-antistop hadroproduction needs a special treatment. At NLO, there is the gq channel as 
an additional production mechanism. At the LHC with a center of mass energy of 7TeV, one 
can expect for a luminosity of 1 fb _1 , 100 to 10 4 events; even 10 5 events are possible, if the stop 
is sufficiently light. At the LHC with a final center of mass energy of 14 TeV, even more events 
are expected to be collected. Hence it is necessary to predict the hadronic cross section with high 
accuracy. 

So far, stop pairs have been searched for at the CDF (SHU and DO experiments I|10H11| at the 
Tevatron using different strategies, for details see Tab. HJ 

Squarks carry colour charge, so it is not surprising that processes involving Quantum Chro- 
modynamics (QCD) obtain large higher-order corrections. For the production of colour-charged 
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superymmetric particles, the NLO corrections have been calculated in Ref. [12], NLL and ap- 
proximated NNLO corrections can be found in the Refs |fT3l - [T6| . It has been found that these 
corrections are quite sizeable. Electroweak NLO corrections to stop-antistop production are dis- 
cussed in Ref. IfTTlHsH . 

The theoretical aspects of tj f*- production up to NLO have been discussed in Ref. J5| and 
of its NLL contributions in Ref. [19|. The hadronic LO and NLO cross section can be evaluated 
numerically using the programme Prospino 11201 . 

In this paper, I calculate and study soft gluon effects to hadronic stop-antistop production in 
the framework of the R-parity- conserving MSSM. I use Sudakov resummation to generate the 
approximated NNLO corrections and include approximated two-loop Coulomb corrections and 
the exact scale dependence. I follow the approach for top-antitop production at the LHC and the 
Tevatron |2lll22| . 

This paper is organised as follows. I review the LO and NLO contributions to the cross section. 
Then I describe the necessary steps to construct the approximated NNLO corrections. Using these 
results I calculate the approximated NNLO cross section and discuss the theoretical uncertainty 
due to scale variation and the error due to the parton density functions (PDFs). I give an example 
how these NNLO contributions reduce the scale uncertainty and improve exclusion limits. 



Ref. 


Process 


Exclusion limit 


Assumptions or comments 


m 


h - 




m h < 100 GeV 


m r o > 50 GeV 


ED 


pp - 


-> hi\ 


130 GeV <m- t < 190 GeV 


Comparison of theor. predictions 
with experimental and observed limits 


EHU 


h - 


> b%f -> hfftPvt 


128 GeV < m h < 135 GeV 




ma 


h - 


> b£+vi 


m h > 180 GeV 
m- h = 100 GeV 


niy > 45 GeV 

75 GeV <niy< 95 GeV 



Table 1: Exclusion limits for stop searches at the Tevatron. 



2 Theoretical Setup 

I focus on the inclusive hadronic cross section of hadroproduction of stop pairs, cr vv ^~ t .~ t * x , which 
is a function of the hadronic center-of-mass energy -^/s, the stop mass m- t , the gluino mass nig, the 
renormalisation scale fl r and the factorisation scale jif. In the standard factorisation approach of 
perturbative QCD, it reads 

s 

^PP/ P p^n*x(s r m^mj^j,^) = J ds Lij(s r s^j) a iHn *(s r m^mj^j r ^) (3) 

i,j=Wg 4m 2 

where the parton luminosities Ly are given as convolutions of the PDFs defined through 

s 

Lij{s,s,vj) = - s J ^f i/p (Vf,^)fj/ V U}, ■ (4) 
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Here, s denotes the partonic center of mass energy and fij, are the factorisation and the renor- 
malisation scale. The partonic cross section is expressed by dimensionless scaling function fj 
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with Ln = In ( ) . The LO scaling functions are given by (5) 
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Formulas for the higher orders of the gg-channel and its threshold expansions can be found in 
Refs IHlQUim, if one takes into account that in the case of stop-antistop production no sum over 
flavours and helicities is needed, 
fit to this function for 
at threshold by 1MH 



(10) 



has been calculated numerically using Prospino ||20| . A 
fit to this function for an easier numerical handling can be found in d. At NLO, f^ 0) is given 
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The constant can be determined from a fit and is approximately given as 0.042 ± 0.001. 

It depends mildly on the squark and gluino masses and on the stop mixing angle. 

The g^-channel is absent at tree level. Its NLO contribution has been extracted from P ro spino. 
This channel is strongly suppressed at threshold. 

The ln /3 terms which appear in the threshold expansions of the NLO scaling functions can 
be resummed systematically to all orders in perturbation theory using the techniques described 
in [24T|28"1) . Logarithmically enhanced terms for the hadronic production of heavy quarks ad- 
mitting an S-wave are also studied in Ref. |29| for arbitrary S!i(3) co i our representations . The 
resummation is performed in Mellin space after introducing moments N with respect to the vari- 
able p = 4m? /s of the physical space: 



a{N,m}) 



1 

J dpp N - l a{s,m}). 
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The resummed cross section is obtained for the individual color structures denoted as I from 
the exponential 



^i(N /m 2 ) 
df jtl (N,ni*) 



g°(m1) ■ exp \G ijA (N + 1)1 + ©(AT 1 log" N) , 
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where all dependence on the renormalisation and factorisation scale ]i r and fir is suppressed and 
the respective Born term is denoted of- j . The exponent Gy x i contains all large Sudakov logarithms 

log^ N and the resummed cross section (fTDl is accurate up to terms which vanish as a power for 
large Mellin-N. To NNLL accuracy G;; i is commonly written as 



G i/;I (N)=logN-4(A)+4 I (A) + ^, 



ii( A ) 



(ii) 
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where A = j6q logNa s / (47r). The exponential exp [G« i(N + 1)] in Eq. ((TTj) is independent from 
the Born cross section (2lll28| . The functions k = 1,2,3, for the octet color structure are 
explicitly given in Ref. |21~| and can be taken over from the case of top-quark hadroproduction, 
the function is given by Eq. (T36)) in the App. IA.2I All g|, < k < 3, depend on a number 
of anomalous dimensions, i.e. the well-known cusp anomalous dimension A q , the functions Dqq 
and Dq controlling soft gluon emission, and the coefficients of the QCD /3-function. The strength 
of soft gluon emission is proportional to the Casimir operator of the SlT(3) co i our representation of 
the produced state. This is identical for tt and £ ; - £*- production. Expressions for An and D„ are 
given in the Refs. |30U31| , and for Dqq in Ref. [32]. At higher orders, they also depend on the 
chosen renormalisation scheme, thus on the dynamical degrees of freedom. 

For my fixed order NNLO calculation, I extracted the a 2 -terms from the right hand side of 
Eq. ((TO)) . At the end, I used Eqs (|30]) - ((35|) given in App. IA.il to convert the Mellin space result 
back to the physical p space. I kept all those terms which are of the order ft 3 \n k ft, < k < 4. 
Eventually, I end up with the following threshold expansion for /i* : 
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is given as = 216 n — 3 2 y, is the unknown 2-loop matching constant, which is 
set to zero in the numerical evaluation and ^qq = 460 — 127T 2 + 72^3 — yf?/, see Ref. [32 J. The 
j6 3 -behaviour of the threshold expansion of the LO cross section comes only from the P-wave of 
the final state £, t* as mentioned in the introduction and does not spoil the factorisation properties 
in the threshold region of the phase space. Note that the formulas given in Ref. [28J can easily 
extended to Mellin transformed cross sections to, which vanish as a power S k with k > 1. These 
two more powers of f> in the ^-channel of £,£*- production lead to an additional 1/N factor 
in Mellin space. Eq. (TO"]) reproduces the known NLO threshold expansion given in Ref. |5l[l9| 
for the qq channel. This is a check that the approach works. Logarithmically enhanced terms 
which are suppressed by an additional 1/N factor appear in the resummation of the (sub)leading 
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log /r (N)/N-terms of the corrections for the structure function Fi and are studied in detail in the 
Refs [33, 34 J. For these reasons I apply the formulas derived for heavy quark hadroproduction. 

The coefficients of the In 4 fi, In 3 fi, and In 2 fi terms depend only on first order anomalous 
dimensions and on the constant ci"* which is related to the NLO constant , see the equation 

above. The linear log fi term depends on as well and on other first order (NLO) contributions, 
but also on second-order anomalous dimensions and non-Coulomb potential contributions |29|. 
In Tab. |2j I show for four examples how these parts contribute to the hadronic NNLO threshold 
corrections. The numbers show that terms which have an NLO origin contribute most and that 
NNLO contributions have a small but sizeable effect. 

I also included the Coulomb corrections up to NNLO. For the singlet case, the Coulomb contri- 
butions are studied in Ref. J35|. Generalisation to other colour structures requires the substitution 
of the corresponding group factors and decomposition of the colour structures of the considered 
process in irreducible colour representations. The last step will not be necessary for stop pair 
production in the ^-annihilation channel. The NLO Coulomb corrections agree with the NLO 
Coulomb corrections for top antitop production (SllMj. In both cases, only the colour octet con- 
tributes to the scaling function at the corresponding leading order in fi. Therefore, I have used 
as an approximation for the NNLO Coulomb contributions for it* -production the same NNLO 
Coulomb contributions as for ^-production |14"ll2Tll3"5ll . Gauge invariance together with Super- 
symmetry support this approximation. Note that the log fi/ /3-term comes from interference of 
the NLO Coulomb-contribution with the NLO threshold logarithms. Tab. |2] shows also the NLO 
and the pure NNLO Coulomb contributions to the NNLO threshold corrections at the hadronic 
level. 

The scale-dependent scaling functions are derived by renormalisation group techniques fol- 
lowing Refs I36U37I : 

4 1 " = ^(wr-c^ 1 -^^ 1 ). <i3 » 

4 21) - (i6^ ( 2 "'4 W) - /T » p S ' - fP » «$') 

+ is? W"' - c ® p " ' - ^ k O • <14) 

^' - (isW ( /r 8 p * ' 8 + 5 /r ' 8 p " r> 8 p ™ + \w ^ p « ' 

+3/44 001 - |wT « <» - \mP » plf) , <i5) 

where (g> denotes the standard Mellin convolution; these are ordinary products in Mellin space 
using Eq. (O. Repeated indices imply summation over admissible partons. However, I restrict 
myself for phenomenological applications to the numerically dominant diagonal parton channels 
at two-loop. Note that the scale dependence is exact at all energies, even away from threshold, 
because the Eqs (fL3l-([T5l) depend on functions which are at least one order lower than they them- 
selves have. The functions P«(x) are called splitting functions and govern the PDF evolution. 
They have the expansion 

p, ( x) = ^p(») (x) + (^) 2 P<;) w + .... (16) 

Explicit expressions for the pjp can be found in Refs I30U38I . 
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Analytical results for fgg 1 ^ and f^' are given in Ref. |23|. For the g^-channel, Eq. ((131 ) simpli- 
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where the functions I/?, L4, and Lg ]23| are defined as 



L 2 = l°g(3). L 4 = U 2 (^)-Li 2 I^E), L„ = log 2 (l-fi)-log 2 (l+fi). 



j±f), L 4 = Li 2 -Li 2 (!±f 
The high energy limit of this scaling function is 
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which agrees with the result given in Ref. 0Q. 

The threshold expansions of the NNLO-scale-dependent scaling functions of the qq channel 
read 
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with d 1 ) = 54na\ 
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Collider 


m~ t 




ln 4 (£) 


ln 3 (/3) 


ln 2 (/3) 


nC 


log(i6) 

D (2) - 
QQ 


rest 


Cnlo 


Cnnlo 


LHC 14TeV 


300 GeV 


74.54 


5.33 


18.80 


31.87 


-2.23 


27.71 
17.90 


12.04 


-9.45 


0.29 


LHC 14TeV 


600 GeV 


2.93 


0.24 


0.81 


1.28 


-0.08 


0.97 
0.63 


0.42 


-0.37 


-0.01 


LHC 7TeV 


300 GeV 


17.5 


1.42 


4.83 


7.66 


-0.47 


5.85 
3,78 


2.54 


-2.21 


-0.06 


Tevatron 


300 GeV 


1.41 


0.16 


0.48 


0.62 


-0.03 


0.34 
0.22 


0.15 


-0.17 


-0.02 



Table 2: Individual hadronic contributions of the log-powers and Coulomb-corrections to the NNLO 

^.(00) 

threshold contributions of the ^-channel in fb. In 4 (j6) has to be understood as fij^zp ' In 4 Q6), anal- 
ogously for the other terms. The linear log-term is decomposed into contributions coming from non- 

(2) 

Coulomb potential terms, from the two loop anomalous dimension D f y q, and, finally, the rest. The 
Coulomb contribution are decomposed into contributions coming from the interference of NLO thresh- 
old logarithms with NLO Coulomb corrections Cnlo an d pure NNLO Coulomb corrections Cnnlo- £ 
denotes the sum over all NNLO threshold contributions. The PDF set used is MSTW 2008 NNLO l39l . 

In Fig. |2j I show the LO, NLO, and NNLO scaling functions. The scaling functions f^°\ fjjl 

f^°\ and depend only on the dimensionless variable tj = ^ — 1, but fQ ' and fffl depend 

also mildly on the masses of the squarks and the gluino and the stop mixing angle (5). At the 
hadronic level, the effect for the NLO + NLL cross section is smaller than 2% [19|, so I neglect 
them. 

As example point, I have chosen the following masses: mj = 300 GeV, rriq = 400 GeV = 
133m- ti , m~ t2 = 480 GeV = 1.6 my , nig = 500 GeV = 1.67 mr, and 9 = n/2, i.e. mj = m- tR 
and m~ ti = m~ t . When varying the stop mass I conserve these mass relations. I restrict myself 
to the lighter stop, but the results also apply to the heavier stop, because the gluon-stop-stop 
interactions entering my process do not distinguish between the left-handed and the right-handed 
stop squarks. 

3 Results 

3.1 Hadronic cross section 

I start with the discussion of the total hadronic cross section, which is obtained by convoluting 
the partonic cross section with the PDFs, see Eq. (O. I keep the gg and the qq channel at all orders 
up to NNLO, and for the scale-dependent terms, only contributions coming from diagonal parton 
channels are considered. Only the NLO contributions of the gq channel are considered, which 
are the leading contributions of this channel. The scale dependence of this channel is given by 
Eq. CLZ). 

I define the NLO and NNLO K factors as 

^NLO 0.NNLO 
K NLO = -^LO~, X NNLO = ^NLO ' ( 23 ) 
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Figure 2: Scaling functions fffi with i = 0,1,2 and j < i. The masses are = 300 GeV, niq = 400 GeV, 
m h = 480 GeV m g = 500 GeV. 
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I use for all orders MSTW2008 NNLO PDFs [39 1, if not otherwise stated. Therefore, the K factors 
account only for the pure higher order corrections of the partonic cross section (convoluted with 
the PDFs) and not for higher order corrections of the PDFs and the strong coupling constant a s . 

In the left column of Fig. [3[ I show the total hadronic cross section for the LHC (7 TeV first 
row, 14 TeV second row) and the Tevatron (third row) as a function of the stop mass. Similar to 
top-antitop |211 l22l and squark-antisquark production [14j, the total cross section shows a strong 
mass dependence. At the LHC (14 TeV), the cross section decreases within the shown stop mass 
range from about 1000 pb to 10~ 2 pb. In the right column of Fig. [H I show the corresponding K 
factors. 

For example, for a stop mass of 300 GeV produced at the LHC with 14 TeV, I have a total cross 
section of 6.57 pb, 9.96 pb, 10.92 pb at LO, NLO, NNLO a pp r ox/ respectively. For a stop mass of 
600 GeV, I find 0.146 pb, 0.216 pb, 0.244 pb. The K-factors are X NLO « 1.5 and Knnlo « 1.1 and 
1.13. For the LHC at 7TeV, one finds similar values for stop masses in the interval 100 GeV < 
m- t < 600 GeV. At the Tevatron, I have Knlo = 1.3 . . . 1.4 and Knnlo ~ 1-2 for stop masses in the 
range 100 GeV < m } < 300 GeV. 

In Tabs [3] - [5[ see App. [BJ values for the total hadronic cross section for different masses, 
PDF sets, scales and colliders are shown. The values for the PDF sets Cteq6.6 [40 1, MSTW 2008 
NNLO 11391 .and CT10 Ell show only small differences, whereas the ABKM09 NNLO (5 flavours) 
PDFs differ in the treatment of the gluon PDF from the other PDF sets. This leads to sizeable 
differences in the total cross sections. 

3.2 Theoretical and Systematic Uncertainty and PDF error 

In this section, I address the following sources for errors: the systematic theoretical error, the scale 
uncertainty, and the PDF error. 

In Tab.[2[ I listed the individual ln^ f> contributions to the total NNLO contributions. Note that 

(2) 

the NNLO matching constants C\- are unknown and set to zero. Compared to the total NNLO 

contributions the In 1 ft term is quite sizeable, this translates to a roughly 3 — 5% contribution to the 
NNLO cross section. To estimate the systematic error coming from the NNLO matching constants 
C\j , I proceed as described in Ref. [22J. I find for the ratio Cnll+CouI / Cexact = 1-10. . . 1.25. This 
ratio translates to an estimate for the relative systematic error coming from the NNLO matching 
constants as 1 — 2.5%. 

The total hadronic LO, NLO, and NNLO cross sections are shown on the left of Fig. 0] as a 
function of the stop mass and for variations of the scale }i with nij/2 < ]i < 2m~ t , where I have 
identified the factorisation scale with the renormalisation scale. The width of the band indicates 
the scale uncertainty, which becomes smaller when going from LO to NLO and NNLO. On the 
right-hand side of Fig. HJ the scale dependence for the example point is shown in more detail. I 
quote as theoretical uncertainty 

mincr(^) < a(nij) < maxcr(ji), (24) 

where the min and max are to be taken over the interval [mj /2,2mj]. At LO and NLO, the 
minimal value is taken at ]i = 2mj and the maximal value at ]i = m~ t /2. However, this is not 
longer true at NNLO. For the theoretical error, one finds 

<r LO = 6.57+2-06 pb/ aNLQ = 9.96t\H pb, <r NNLO = 10.90±g;?| pb . (25) 

As one can see, there is a strong scale dependence at LO, becoming weaker at NLO, and is 
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Figure 3: Total hadronic cross section at LO, NLO, and NNLO a pp r0 x at the LHC 7TeV(first row) and 
14TeV(second row) and the Tevatron (1.96 TeV, third row). The right column shows the corresponding K 
factors. The PDF set used is MSTW2008 NNLO (39). 

flattend out at NNLO within the considered range. This flattening gives a hint that the approach 
is reliable. 

Using renormalisation group techniques, one recovers the full dependence on the renormal- 
isation scale \i r and factorisation scale ]df. I have done this for the example point for the NLO 
and the NNLO cross section, see Fig. |5l I define the theoretical uncertainty coming from an 
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Figure 4: Left hand side: Theoretical uncertainty of the total hadronic cross section at the LHC (14 TeV) at 
LO (upper figure, blue band), NLO (central figure, green band), NNLOapprox (lower figure, purple line). 
At NNLO a pp r ox, the theoretical uncertainty has shrunk to a small band. Right hand side: Scale dependence 
of the total hadronic cross section for the example point m% — 300 GeV, niq = 400 GeV, nij 2 = 480 GeV, 
nig = 500 GeV. The vertical bars indicate the total scale variation in the range [mj/2,2mj]. 



(26) 



independent variation of }i r and fit in the standard range }i r ,}if € [rtif /2,2mj] as 

mincer, fif) < cr(mi) < maxc r (^ r , }if). 

The contour lines of the total cross section for the example point with an independent variation 
of }i r and jif are shown in Fig. |5j Note that the range of the axes is from log 2 (^ r j/wj) = — 1 
to log 2 (/i r f/tn~ t ) = 1. The scale variation with fixed scales y, r = ]it proceeds along the diagonal 
from the lower left to the upper right corner of the figure. The gradient of the NLO contour 
lines lies approximately in the ]i r = fit direction, meaning that the theoretical error from the 
definition in Eq. (|26l is the same as if one sets ]i r = jif. For the NNLO case one observes the 
opposite situation: the contour lines are nearly parallel to the diagonal ji r = jif. I obtain a larger 
uncertainty in that case: 



PiMLO = 9 -%-i22 P b ' ^NNLO 



10 Qf) +1 - 05 



pb. 



(27) 



Another source of error to discuss is the PDF error. I calculated the PDF uncertainty according 
to Ref. HO for the two PDF sets CT10 and MSTW2008 NNLO (90% C.L.). In both cases, the 
uncertainty increases with higher stop masses due to large uncertainties of the gluon PDF in 
high x-ranges. For CT10, I find as relative errors ~ 3% for m- t = 100 GeV and ~ 18% for mj = 
1000 GeV, and for MSTW2008 NNLO, the relative errors are w 3% for m- t = 100 GeV and w 10% 
for m- t = 1000 GeV. The relative error of the MSTW2008 NNLO PDF set is smaller for large stop 
masses compared to the CT10 PDFs. 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

log2(/i_r/m_stop) log2(/i_r/m_stop) 

Figure 5: Contour lines of the total hadronic NLO (left) and NNLO (right) cross section from the inde- 
pendent variation of the renormalisation and factorisation scale \i r and \i f for LHC, 14 TeV, with PDF set 
MSTW2008 NNLO (39) for the example point with m~ t = 300 GeV. The dot in the middle of the figure 
indicates the cross section for \i r = \ij — m~ t , and the range corresponds to }i r 6 [m- t /2,2nif\. 
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Combining theoretical uncertainty and PDF error one obtains 

(7 NNL o = 10.90 pb loi8P b ( scale ) ±a55P b (MSTW2008 NNLO) (28) 
tfNNLO = 10.86 pb ig;gpb (scale) ±°;gpb (CT10). (29) 

3.3 Mass exclusion limits 

The approximated NNLO contributions enlarge the stop-antistop production cross section by 
10 — 20%, depending on the hadron collider, its center of mass energy, and the stop mass. This 
could be converted into larger exclusion limits for the mass of the stop squark. A stop with a mass 
my = 120 GeV has a NLO production cross section of 5.05 pb, the same cross section corresponds 
at NNLO to a stop with a mass of 123.5 GeV. At higher stop masses, the exclusion limit is even 
further enhanced: 0.164 pb corresponds at NLO to a stop with a mass of 210 GeV, but at NNLO 
to a mass of 215 GeV. At the LHC (14 TeV) one would have a similar situation. An NLO cross 
section of 750 pb corresponds to m~ t = 120 GeV, but to m~ t = 121.5 GeV at NNLO. And if the stops 
are heavier, a cross section of 10 pb is related to mi = 300 GeV at NLO, but to mr = 305.5 GeV at 
NNLO. 

4 Conclusion and Summary 

In this paper, I computed the NNLO threshold contributions including Coulomb corrections for 
stop-antistop production at hadron colliders. 

• I presented analytical formulas for the threshold expansion of the NNLO scaling function 
using resummation techniques for the scale independent scaling function and RGE tech- 
niques for the scale-dependent scaling functions. 

• After convolution with suitable PDF sets, the NNLO corrections are found to be about 20% 
for the Tevatron and 10 — 20% for the LHC compared to the hadronic NLO cross section. The 
PDF sets Cteq6.6 (H, MSTW 2008 NNLO (H, and CT10 EU show only small differences 
in the total cross section, whereas the values obtained with the PDF set ABKM09 NNLO (5 
flavours) [42 J differ by 10 — 35% due to differences in the gluon PDF. 

• I calculated the exact scale dependence and I found a remarkable stabilisation of the cross 
section under scale variation. For my example point, the theoretical error is reduced from 
12% at NLO to better than 2% at NNLO. 

• I discussed three types of errors: systematic theoretical errors, uncertainties due to scale 
variation and PDF errors. The systematic error was estimated to be about 3 — 6%, the scale 
uncertainty to be about 2% or better, and the PDF error 3 — 18% depending on the stop mass 
and the PDF set used. 

• Finally, I demonstrated how NNLO cross sections could enlarge exclusion limits. The im- 
provement of the lower exclusion limit were about a few GeV. 
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A Useful Formulas 

A.l Mellin Transformations 



o 

1 

J dp^-yin 3 /} 



dp^-^ln 4 /* 



dpp N -^ 2 ]n^ 



dpp N -^ 2 ]n 2 ^ 



8 + 61n2 + 31nN 



(30) 



32 



3/r 2 + 16 - 641n 2 + 241n 2 2 + 8( - 4 + 31n 2) In N + 61n 2 N 
1 



(31) 



3v^r 



64 



256 



32 In 2 - 64 In 2 2 + 16 In 3 2 - 8n 2 + 6/r 2 In 2 + 28£ 3 
- (16 + 3/r 2 - 64 In 2 + 24 In 2 2) In N 

-2 ( - 8 + 6 In 2) In 2 N + 2 In 3 n] -J^ (l + (1/N) ) (32) 



21tt 4 + 12/r 2 (8 - 32 In 2 + 12 In 2 2) 
+4 (192 In 2 2 - 256 In 3 2 + 48 In 4 2 - 448£ 3 + 336£ 3 In 2) 
+24( - 8/r 2 + 6/r 2 In 2 + 32 In 2 - 64 In 2 2 + 16 In 3 2 + 28£ 3 ) In N 
+12(16 + 3/r 2 - 64 In 2 + 24 In 2 2) In 2 N 



16(-8 + 61n2)ln 3 N + 121n 4 N -J^ (l + O (1/N) ) (33) 

(34) 



ilnNl(l + 0(l/N)) 



1 r 
24 



7T 2 - 121nN + 61n 2 N 



W (l + 0(1/N)) 



(35) 



with lnN = lnN + 7 e . A detailed overview of Mellin transformations is given in Ref. 11431 
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A.2 Resummation 

For the sake of clarity, I specify the function g^ 8 I have used in my resummation. The func- 
tion g°. g differs a little bit from the function given in Ref . [21 J due to Mellin transformation of 
functions of the type /5 3 log' c (/3) instead of /31og' c (/3), < k < 4, as it is in the case of ££-production. 



£qq,8 



1 + Cis 



C F ( - 64 + 8j e 2 + 4tt 2 - 321n 2 + 81n 2 2) + C A (-8 + 4j e - 41n2) + 4C, 
-2C F (16-87 e )ln2 



(i) 
<?<? 



C F C 



(i) 
m 



3328 



+ 327 e 2 + 16tt 2 + 512 In 2 - 647 e In 2 - 256 In 2 2 



^ / 30592 448, 224 



27 



80 2 32 , 56 2 14272, „ 

— ^3 " ^le - ~^Te - -^-Te ~ ~ 71 " ^^ln2 



56 

3 



27 



160 
~9~ 



I_l 7r i n 2 + f 7 2 In 2 + 16/r 2 In 2 + ^ In 2 2 - °-^ le In 2 2 1 - 



32 

T 



832 



ln J 2 



-C 



335872 14336 



4352 



3 £ 3 - 512 7 ^ + 32 7 * + ^- tt z + 327i 2 7 2 + 24/r 1 



+ 



27 
57344 



In 2 - 5376^3 In 2 + 10247 e In 2 - I287 3 In 2 - 1024/r 2 In 2 



-647T z 7 e In 2 - 15872 In 2 2 + I927 2 In 2 2 + 6087T 2 ln z 2 + 9216 In" 3 2 



.2 i„2 , 



-1287 e In 3 2 -2560 In 4 2 



f-^ + 167, + 32 In 2 



/512 88 8 2 4 2 464, „ 16 , „ 64, 2 \ 
" C ^ n / ( -9- - "9-^ - - 3^ 2 -— In2 + - 7e ln2 + -ln 2 2j 



+C A C F 



250432 7840 „ 5296 



81 



+ 



40 



464 , 5188 



27 



-71 



-16/r z 7 e - ^/r 2 7 2 - ^ + 2 „ 

464 2 , „ 728 2 , ^ 16 2 , n 33088, 2 „ 464 , 2 „ 
-— 7 2 In 2 — 7i 2 In 2 + — n 2 j e In 2 — In 2 2 + —y e In 2 



132256, „ „ nT , „ 80 , „ 
In 2 - 112£ 3 In 2 + y 7e In 2 



64 2l 2 „ 12064, 3 „ 
+ y7T 2 ln 2 2+ — — In 3 2 



+C 



2816 172 



134 2 4 -> 



3608 



H 7e H 7» H 7T 7T 7e + 

9 9 ' 3 £ 9 3 ' 9 

I ^2 m 2 - ^ In 2 2 - + I6C3 In 2 + 87^3 



1 n 136 1 „ 

ln2- — 7,ln2 



-C 



(2) 



(36) 



7 e denotes the Euler constant and £3 is the Riemannian Zeta function evaluated at x = 3. 
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B Numerical Results 





o-(LO)[pb] 


^(NLO)[pb] 


cr(NNI 


X>)[pb] 


[GeV] 


x = A x = 1 x = 2 


x = i_ X = 1 x = 2 


i 

x = J x 


= 1 x = 2 



Cteq6.6 



120 


579.745 


461.516 


372.340 


842.538 


741.438 


652.012 


762.236 


786.927 


785.150 


180 


97.170 


75.718 


60.001 


133.421 


118.340 


103.980 


124.929 


127.242 


126.175 


240 


25.299 


19.463 


15.254 


33.578 


29.929 


26.279 


32.176 


32.490 


32.086 


300 


8.466 


6.455 


5.020 


10.976 


9.817 


8.615 


10.702 


10.739 


10.573 


360 


3.335 


2.526 


1.953 


4.250 


3.810 


3.342 


4.201 


4.196 


4.121 


480 


0.710 


0.532 


0.408 


0.883 


0.794 


0.696 


0.892 


0.884 


0.865 


600 


0.198 


0.147 


0.112 


0.242 


0.218 


0.191 


0.249 


0.245 


0.239 


MSTW 2008 NNLO 


120 


587.154 


468.209 


378.090 


849.643 


749.028 


659.621 


769.079 


793.744 


791.998 


180 


99.033 


77.177 


61.136 


135.263 


120.146 


105.649 


126.755 


129.046 


127.921 


240 


25.825 


19.848 


15.538 


34.083 


30.413 


26.712 


32.698 


32.996 


32.566 


300 


8.632 


6.571 


5.101 


11.127 


9.960 


8.740 


10.867 


10.896 


10.718 


360 


3.390 


2.562 


1.977 


4.295 


3.854 


3.379 


4.255 


4.245 


4.165 


480 


0.714 


0.534 


0.408 


0.883 


0.795 


0.696 


0.895 


0.886 


0.866 


600 


0.197 


0.146 


0.111 


0.239 


0.216 


0.189 


0.247 


0.243 


0.236 


CT10 


120 


584.754 


465.703 


375.794 


849.755 


747.907 


657.825 


768.385 


793.474 


791.811 


180 


98.265 


76.560 


60.641 


134.841 


119.609 


105.071 


126.243 


128.588 


127.483 


240 


25.605 


19.688 


15.415 


33.955 


30.263 


26.555 


32.541 


32.854 


32.425 


300 


8.568 


6.528 


5.070 


11.098 


9.923 


8.699 


10.824 


10.858 


10.679 


360 


3.374 


2.553 


1.970 


4.294 


3.848 


3.371 


4.248 


4.239 


4.158 


480 


0.717 


0.537 


0.410 


0.890 


0.800 


0.700 


0.900 


0.891 


0.870 


600 


0.199 


0.148 


0.112 


0.244 


0.219 


0.191 


0.250 


0.246 


0.240 


ABKM 09 NNLO (5 flv) 


120 


539.054 


428.594 


345.556 


764.660 


677.764 


598.338 


700.513 


720.046 


718.084 


180 


87.613 


68.261 


54.112 


117.982 


105.305 


92.854 


111.788 


113.439 


112.487 


240 


22.132 


17.046 


13.380 


28.906 


25.912 


22.833 


28.027 


28.207 


27.868 


300 


7.190 


5.498 


4.288 


9.196 


8.270 


7.284 


9.075 


9.079 


8.947 


360 


2.752 


2.093 


1.626 


3.464 


3.123 


2.751 


3.468 


3.454 


3.397 


480 


0.554 


0.419 


0.323 


0.682 


0.617 


0.544 


0.698 


0.690 


0.677 


600 


0.146 


0.110 


0.085 


0.177 


0.161 


0.142 


0.185 


0.182 


0.178 



Table 3: Numerical values for the stop pair-production cross section at LHC with -y/s = 14TeV and the 
PDF sets Cteq6.6 (40), MSTW 2008 NNLO (39), ABKM09 NNLO (5 flavours) 021, and CT10 [41]. The QCD 
predictions are given at LO, NLO, and NNLO accuracy and for different stop masses and scales x = }i/mj. 
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cr(LO)[pb] 


tj(NLO)[pb] 


t7(NNLO)[pb] 


[GeV] 


x = 1 x = 1 x = 2 


x = 2 X = 1 x = 2 


x = 2 x = 1 x = 2 



Cteq6.6 



120 


134.391 


101.197 


77.851 


184.136 


160.840 


138.709 


173.624 


176.078 


173.026 


180 


18.038 


13.341 


10.106 


23.542 


20.725 


17.855 


23.084 


23.068 


22.511 


240 


3.881 


2.839 


2.130 


4.925 


4.354 


3.748 


4.957 


4.907 


4.765 


300 


1.091 


0.792 


0.590 


1.359 


1.205 


1.036 


1.395 


1.372 


1.327 


360 


0.364 


0.263 


0.195 


0.448 


0.398 


0.342 


0.468 


0.457 


0.441 


420 


0.137 


0.098 


0.072 


0.168 


0.149 


0.128 


0.177 


0.172 


0.165 


MSTW 2008 NNLO 


120 


137.257 


103.302 


79.392 


186.837 


163.495 


141.098 


176.403 


178.825 


175.655 


180 


18.377 


13.561 


10.249 


23.825 


20.998 


18.086 


23.416 


23.378 


22.789 


240 


3.918 


2.857 


2.137 


4.941 


4.371 


3.759 


4.991 


4.933 


4.782 


300 


1.087 


0.786 


0.584 


1.346 


1.194 


1.025 


1.388 


1.362 


1.315 


360 


0.357 


0.257 


0.190 


0.437 


0.388 


0.333 


0.458 


0.447 


0.430 


420 


0.132 


0.094 


0.069 


0.160 


0.142 


0.122 


0.170 


0.165 


0.159 


CT10 


120 


136.074 


102.422 


78.752 


186.261 


162.710 


140.298 


175.639 


178.121 


175.005 


180 


18.262 


13.496 


10.211 


23.805 


20.956 


18.042 


23.354 


23.332 


22.753 


240 


3.922 


2.866 


2.146 


4.971 


4.394 


3.777 


5.008 


4.955 


4.805 


300 


1.100 


0.797 


0.592 


1.368 


1.212 


1.040 


1.406 


1.381 


1.333 


360 


0.366 


0.263 


0.195 


0.450 


0.399 


0.342 


0.470 


0.459 


0.441 


420 


0.137 


0.098 


0.072 


0.167 


0.148 


0.127 


0.177 


0.172 


0.165 


ABKM 09 NNLO (5 flv) 


120 


117.501 


88.531 


68.185 


157.745 


138.823 


120.224 


151.054 


152.485 


149.896 


180 


14.825 


11.007 


8.373 


19.061 


16.894 


14.626 


18.991 


18.902 


18.474 


240 


3.004 


2.215 


1.674 


3.765 


3.353 


2.904 


3.856 


3.804 


3.704 


300 


0.797 


0.586 


0.441 


0.981 


0.877 


0.761 


1.027 


1.006 


0.978 


360 


0.252 


0.185 


0.139 


0.306 


0.274 


0.238 


0.326 


0.318 


0.308 


420 


0.090 


0.066 


0.049 


0.108 


0.097 


0.084 


0.117 


0.113 


0.110 



Table 4: Numerical values for the stop pair-production cross section at LHC with y/s = 7 TeV and the same 
PDF sets as in Tab. [3] The QCD predictions are given at LO, NLO, and NNLO accuracy and for different 
stop masses and scales X = ji/mj. 
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m~ t 


a 


(LO)[pb] 




a 


!NLO)[pb] 


o-Qs 


INLO)[] 


pb] 


[GeV] 


x- i 

X — 2 


x = 1 


x = 2 


x- i 

X — 2 


x = 1 


x = 2 


x - i 

X — 2 


x = 1 


x = 2 


Cteq6.6 


120 


5.525 


3.850 


2.779 


5.782 


5.246 


4.509 


6.453 


6.144 


5.780 


150 


1.570 


1.095 


0.790 


1.572 


1.448 


1.254 


1.793 


1.695 


1.589 


180 


0.525 


0.366 


0.264 


0.519 


0.480 


0.417 


0.597 


0.562 


0.524 


210 


0.194 


0.135 


0.097 


0.192 


0.178 


0.154 


0.222 


0.208 


0.193 


240 


0.076 


0.053 


0.038 


0.076 


0.070 


0.061 


0.088 


0.082 


0.076 


270 


0.031 


0.021 


0.015 


0.032 


0.029 


0.025 


0.036 


0.034 


0.031 


300 


0.013 


0.009 


0.006 


0.013 


0.012 


0.010 


0.015 


0.014 


0.013 


M 


STW 2008 NNLO 


120 


5.336 


3.712 


2.677 


5.562 


5.047 


4.336 


6.242 


5.930 


5.568 


150 


1.492 


1.040 


0.751 


1.482 


1.368 


1.186 


1.701 


1.605 


1.503 


180 


0.493 


0.344 


0.248 


0.482 


0.448 


0.390 


0.558 


0.524 


0.489 


210 


0.180 


0.126 


0.090 


0.177 


0.164 


0.143 


0.205 


0.191 


0.178 


240 


0.070 


0.049 


0.035 


0.069 


0.064 


0.056 


0.080 


0.075 


0.069 


270 


0.028 


0.020 


0.014 


0.029 


0.026 


0.023 


0.033 


0.031 


0.028 


300 


0.012 


0.008 


0.006 


0.012 


0.011 


0.009 


0.014 


0.013 


0.012 


CT10 


120 


5.546 


3.861 


2.784 


5.793 


5.256 


4.516 


6.479 


6.164 


5.795 


150 


1.570 


1.094 


0.789 


1.567 


1.444 


1.251 


1.792 


1.693 


1.586 


180 


0.524 


0.365 


0.263 


0.515 


0.478 


0.415 


0.595 


0.559 


0.521 


210 


0.193 


0.134 


0.096 


0.190 


0.176 


0.153 


0.220 


0.206 


0.191 


240 


0.075 


0.052 


0.037 


0.075 


0.070 


0.060 


0.087 


0.081 


0.075 


270 


0.031 


0.021 


0.015 


0.031 


0.029 


0.025 


0.036 


0.033 


0.031 


300 


0.013 


0.009 


0.006 


0.013 


0.012 


0.010 


0.015 


0.014 


0.013 



ABKM 09 NNLO (5 fly) 



120 


4.478 


3.207 


2.371 


4.374 


4.095 


3.617 


5.023 


4.773 


4.532 


150 


1.312 


0.942 


0.698 


1.222 


1.165 


1.040 


1.423 


1.348 


1.278 


180 


0.453 


0.325 


0.240 


0.419 


0.401 


0.359 


0.488 


0.461 


0.435 


210 


0.172 


0.123 


0.090 


0.161 


0.153 


0.137 


0.186 


0.175 


0.165 


240 


0.069 


0.049 


0.036 


0.066 


0.062 


0.055 


0.076 


0.071 


0.067 


270 


0.028 


0.020 


0.015 


0.028 


0.026 


0.023 


0.032 


0.030 


0.028 


300 


0.012 


0.008 


0.006 


0.012 


0.011 


0.010 


0.014 


0.013 


0.012 



Table 5: Numerical values for the stop pair-production cross section at the Tevatron with \/s = 1.96 TeV 
and the same PDF sets as in Tab. [3] The QCD predictions are given at LO, NLO, and NNLO accuracy and 
for different stop masses and scales X = ji/utj. 
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